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ABSTRACT 

It is argued that depolarization of solar radio bursts requires reflection off 
boundary layers no thicker than about a wavelength (a few meters at most) 
between regions with large density ratios. The implied inhomogeneities suggest 
that the corona is much more highly and sharply structured than can be resolved 
from observations at other wavelengths. A simplified version of magnetoionic 
theory is used to derive a depolarization coefficient: the effect of the magnetic 
field is ignored in treating the dispersion, but taken into account in treating the 
(circular) polarization. Plots of the depolarization coefficient are used to infer 
conditions under which effective depolarization occurs, and it is concluded that 
favorable conditions require total internal reflection. For type I sources away from 
the central meridian, effective depolarization requires reflection off an overdense 
structure with a density ratio £ ~ 10. For type III bursts, a density ratio £ ~ 2 
suffices, with at least two reflections off walls of ducts at w 20° to the radial. 

Subject headings: solar corona, solar radio bursts, polarization, inhomogeneities 



1. Introduction 

There is compelling evidence for a depolarizing agent operating on solar radio emission 
as it propagates through the solar corona. The evidence is from depolarization of funda- 
mental (F) plasma emission, which simple theory predicts should be 100% in the sense of 
the o mode (Melrose 1980, 1985, 1991). This prediction is supported by observations of 
type I emission, which is close to 100% polarized in the sense of the o mode for sources 
near the central meridian (Elgar0y 1977). However, type I sources closer to the limb have 
an intermediate polarization that decreases as the source approaches the limb (Zlobec 1975; 
Wentzel, Zlobec and Messerotti 1986). A model for the directivity of type I emission invokes 
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reflection off overdense fibers (Bougeret and Steinberg 1977), for sources away from the cen- 
tral meridian, and depolarization is thought to be associated with this reflection (Wentzel, 
Zlobec and Messerotti 1986). There is also evidence for intermediate polarization in type IV 
emission (Chernov & Zlobec 1995), indicating that this effect is not unique to type I sources. 
Depolarization is required for all F emission in type II and type III bursts. Such F emission 
can be moderately to relatively highly polarized (~ 70%) in the sense of the o mode (Dulk, 
Suzuki & Sheridan 1984; Suzuki and Dulk 1985). However, it is never 100% polarized. If 
depolarization is due to reflection, this implies that type II and type III F emission always 
experiences at least one depolarizing reflection before escaping from the corona. There is 
strong evidence for ducting of both F and H (harmonic) type III emission (Duncan 1979), 
and one suggestion is that the depolarization results from reflection off the walls of a duct 
(Hayes 1985b). 

In this paper it is assumed that depolarization of solar radio bursts is due to reflection 
off field-aligned, overdense structures with sharp boundaries, and the implications of this 
assumption are explored. The theory of reflection and transmission of magnetoionic waves 
at a sharp boundary exhibits a rich and varied range of features (Hayes 1985a,b). This 
makes it difficult to identify the general features needed in formulating a semiquantitative 
theory of depolarization. A relatively simple analytic model for the depolarization is needed. 
Such a model is proposed here: it corresponds to a limit of magnetoionic theory in which 
the magnetic field is assumed to be arbitrarily small. This is a valid approximation in 
cases where the polarization of the natural modes is approximately circular, and when the 
difference between the refractive indices of the o and x modes is small compared with their 
mean. The most serious limitation of this approximation is near the plasma frequency, 
where there is a range of frequencies where the o mode can propagate and the x mode 
is either evanescent or is becomes the z mode. (This limitation is no important provided 
that the reflecting boundary is well above the source region, which is the case in the models 
discussed here.) For incident radiation that is 100% circularly polarized, this model is used to 
calculate the reflection coefficients for radiation with the same handedness and the opposite 
handedness. The net degree of circular polarization of the reflected radiation is identified as 
the depolarization coefficient. The reflection is assumed to be off a sharp boundary over an 
overdense region, with a density ratio £, aligned along magnetic field lines, with the incident 
ray at an angle 6 to the field line. 

The requirement for reflection to lead to depolarization is discussed in section 2: existing 
arguments that the boundary can be no more than about a wavelength (a few meters) thick 
(Budden 1961; Hayes 1985a) are supported by an argument based on mode-coupling theory. 
The simplified form of magnetoionic theory is introduced in section 3 and used to calculate 
a depolarization coefficient for reflection off a sharp boundary. In section 4 it is argued 
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that effective depolarization requires total internal reffection, and favorable conditions for 
depolarization of type I and type III bursts are identified. The results are summarized and 
discussed in section 5. 



2. Reflection versus refraction 

The difference between refraction through a large angle in a boundary layer and re- 
flection at a sharp boundary is discussed in this section. When refraction applies, wave 
propagation in an inhomogeneous medium may be described using mode-coupling theory, 
with each mode described by a WKB solution. Reflection applies when the density gradient 
is so steep that mode coupling theory is invalid. 

2.1. Qualitative discussion of mode coupling 

If any inhomogeneity is sufficiently weak, radiation initially in a single magnetoionic 
mode remains in that mode. A small leakage into the other mode is described by the theory 
of mode coupling. A standard treatment of mode coupling can be summarized as follows 
(Budden 1961; Melrose 1974a,b). Suppose the plasma density (or some other plasma param- 
eter) increases smoothly along a specific direction, identified as the z-axis. The equations 
of the magnetoionic theory are Fourier transformed in time and in x and y, and reduced to 
four coupled linear differential equations in z for E x , E y , B x , B y . In the homogeneous case, 
on also Fourier tranforming in z, the condition for a solution to exist reduces to a quartic 
equation (the 'Booker quartic' equation) for k z . In simple cases, the four solutions may be 
interpreted as four different modes and labeled as of, o|, xf and x j, where the arrows refer 
to the direction of propagation relative to the z-axis. Including the inhomogeneity leads to 
coupling between these four modes. Coupling becomes strong near coupling points, where 
the values of k z for two modes are equal. Leakage between one mode and the other is due 
to mode coupling between o| and xf (or o|, and x|). A reflection point for the o mode 
or the x mode refers to a point along a curved ray path where k z passes though zero and 
a refracting ray reverses its sense of propagation along the direction of the gradient. At a 
reflection point, coupling between o| and o| (or xf and x|) becomes strong. For the theory 
of mode coupling to be valid, the coupling points must be well separated, requiring that the 
coupling between say o j and x j near the reflection point for the o mode or the x mode must 
be weak for the theory to apply. 

In contrast, at a sharp boundary, an incident wave in one mode produces reflected (and, 
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if allowed, transmitted) waves in both modes. Mode coupling theory cannot be valid when 
reflection leads to two reflected modes, which is the case of interest in explaining depolariza- 
tion. The transition from reflection-like refraction to true reflection cannot be treated within 
the framework of mode-coupling theory, and requires a full-wave theory (Budden 1961; Hayes 
1985a). 

2.2. Separation of reflection points 

When the density gradient is sharp, reflection and transmission are treated by applying 
the electromagnetic boundary conditions at the interface. This requires that the reflection 
points, in the sense used in mode-coupling theory, be effectively coincident. This qualitative 
condition can be converted into a semi-quantitative condition for reflection as follows. 

Given an analytic model for the boundary layer, geometric optics allows one to determine 
the location of the reflection points for the o and x modes inside the boundary layer. In the 
vicinity of a reflection point, the full wave solution may be approximated by an Airy function 
(Budden 1961), which is an oscillating function on one side and an exponentially decaying 
function on the other side of the reflection point. At the reflection point for the x mode, 
which is on the lower-density side of the reflection point for the o mode, the x mode becomes 
evanescent. For reflection to apply in this case, the skin depth of the x mode beyond 
its reflection point must be large in comparison with the distance between the reflection 
points for the x mode and the o mode. Only then can radiation tunnel from one reflection 
point to the other before decaying significantly. The way this skin depth is estimated is 
outlined in Appendix A, cf. (A.l). Let the density gradient be characterized by a distance 
L = l/|gradn e |. The requirement for reflection to apply, when the cyclotron frequency is 
much smaller than the plasma frequency, Q e <C u> p , then becomes 

L<(c/cu)[(l-r 2 )(l + r 2 cos 2 0)] 1/2 , (1) 

where is the angle between the incident ray direction and the plane defined by the density 
gradient and B. The parameter r may be interpreted as the sine of the angle of incidence 
times the initial refractive index. It follows that the plasma frequency much change over a 
distance of order the free-space wavelength in order for reflection rather than refraction to 
apply (Budden 1961; Hayes 1985a). 



- 5 - 



3. Depolarization due to reflection 



Depolarization due to reflection off a sharp boundary is due to incident radiation in 
one mode resulting in reflected radiation in both modes. A simplified model is used here to 
derive a formula for the depolarization on reflection. 



The approximate treatment of depolarization due to reflection is based on the following 
assumptions. 

1. The magnetic field is neglected in the refractive indices for the magnetoionic modes 
but is taken into account in their polarizations, which are assumed circular, denoted 



2. The boundary layer is parallel to the magnetic field lines. 

3. The incident wave is in the plane defined by the magnetic field lines and the normal 
to the boundary layer. 

The following parameters need to be specified. The ratio of the wave frequency to the plasma 
frequency on the low-density side is specified by X = uj 2 /uj 2 , and the corresponding ratio 
on the high-density side by X' = uj' 2 /uj 2 = £X. The angle of incidence also needs to be 
specified, and it is convenient to write this as 9 mc = n/2 — 9, where 9 is the angle between 
the incident ray and the magnetic field. Let the incident ray be in a plane at an azimuthal 
angle 0. In general, there is a transmitted ray, at angles 9', 0' and a reflected ray at 9", 0". 
These satisfy Snell's law in the form 



Assumption 3 corresponds to = 0, with the incident, reflected and transmitted waves all in 
the plane containing the direction of the gradient and the magnetic field. Then the reflected 
ray is deflected through an angle #def = 20. 

In this model, there are only two different reflection coefficients: r RR = r LL and r RL = 
r LR , where the first (second) subscript indicates the polarization of the incident (reflected) 
mode. Assuming incident i?-polarization, the depolarization coefficient is 



3.1. Simplified model for reflection coefficients 



R,L. 



(1 -X') 1/2 cos9' = (1 -X) 1/2 cos9, 0' = 0, 9" = 9, 0" = + vr. 



(2) 



rRR - r RL 



(3) 



PR = 



trr + r RL 
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Calculation of the Fresnel coefficients in this model involves modifying the standard text- 
book derivation by separating the incident (unprimed), transmitted (primed) and reflected 
(double primed) waves into circularly polarized components, E Rj l, E' rl , E rl . For inci- 
dent R polarization one sets El = 0, and solved for the Fresnel coefficients E RL /E R , with 
r RR = \E R /E R \ 2 t R l = \E'[/E R \ 2 . The total reflection coefficient is r tot = r RR + r RL , and 
the transmission coefficient is 1 — r tot . 

Three different cases need to be distinguished in writing down the Fresnel coefficients: 

A. For X < X' < 1 — (1 — X) cos 2 0, there are both transmitted and reflected modes. 

B. For 1 — (1 — X) cos 2 9 < X' < 1, there is total internal reflection and no transmitted 

wave. As in a standard derivation of the Fresnel coefficients, one can use the result 
in case A by formally allowing (in this case) cos0' > 1, with sin0' = (1 — cos 2 ^') 1 / 2 
imaginary. 

C. For X' > 1 the refractive index in the denser medium is imaginary, and there is no 

transmitted wave. In this case one has cos0' imaginary, with sin0' = (1— cos 2 0') 1 / 2 > 1. 



It is convenient to introduce parameters 

a = sin 0'/ sin 0, c = cos0/ cos0', (4) 

which are real in case A and one of which is imaginary in cases B and C, respectively. One 
may rationalize the expressions for the Fresnel coefficients such that a, c appear only either 
squared or in the combination ac. Then in both cases B and C ac is imaginary, and these 
two cases can be treated together. Before rationalization, the Fresnel coefficients are 

K = c(l-a 2 ) El = -q(l - c 2 ) 

E R (a + c)(l + ac)' E R (a + c)(l + ac)' 



3.2. Reflection and transmission 

The reflection coefficients may be calculated directly for case A using (5). One finds 

cos 29 + cos 29' nn . X' - X + (1 - X) cos 29 
PR = 1 + cos 20 cos 20" C ° S 29 = T=X> • (6) 

The total reflection coefficient is 

(1 + cos 20 cos 20') (cos 20 - cos20') 2 
rtot ~ 2 sin 2 (0 + 0') (sin 20 + sin 20') 2 ' ( ' 
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The angle 9' increases with increasing frequency, from 9' = at (uj 2 — uj 12 ) sin 2 9 = uj' 2 — uj 2 
to 9' — > 9 for uj — > oo. 

For X' 1 most of the energy is transmitted (Hayes 1985a), so that this case is of little 
interest here. The degree of depolarization written down by Wentzel, Zlobec and Messerotti 
(1986) is reproduced by the high-frequency limit, X, X' — > 0, of (6), when one has 9' = 9 
and 29 is identified as the angle of deflection. A different derivation of this limiting case is 
given in Appendix B, and this alternative derivation shows that the result applies when the 
incident ray is not in the plane defined by the magnetic field and the normal to the boundary. 



3.3. Total internal reflection 

In cases B and C there is total internal reflection, so that one has r tot = I. After 
rationalization one finds 

(c 2 + a 2 )(l + c 2 a 2 )-Ac 2 a 2 
PR ~ (c 2 - a 2 )(l - c 2 a 2 ) ' ( ) 

with either a 2 < (case B) or c 2 < (case C). In case C, c 2 is negative, and it is convenient 
to write c 2 = -1/x, with x = (1 - X)/(X' - 1) = (uj 2 - uj 2 p )(uj' 2 - uj 2 ). One then has 
a 2 sin 2 9 = 1 + x cos 2 9 > 0. Thus in case C, (8) simplifies to 

x — (1 + 3x) sin 2 9 + 2x sin 4 9 
x + (l-x) sin 2 9 

One finds that (9) also applies in case B, where a 2 and x < — 1/ cos 2 9 are negative. 



PR = —T, 3 - 2 n ' ( 9 ) 



The two expressions (6), for case A, and (9), for cases B and C, cover the entire range of 
angles and frequencies, uj > uj p . The boundary between their regions of validity is 9' = in 
(6) and x — — 1/ cos 2 9 in (9). Contour plots of pn as a function of uj/uj p and 9 have similar 
forms for different r, provided the axes are scaled appropriately. An approximation to (9) 
that exhibits this scaling is for x 1 and sin 2 9 <^ 1 when one has 

1 - 3£(u6 1 W) 2 , s 

1 + £{u9/u p ) 2 

which applies for 1 < oj/jj p < £ 1/2 , < 1. 



3.4. Depolarization coefficient 



Numerical results for the depolarization coefficient are shown as a contour plot in fig- 
ure 1. The depolarization coefficient is close to unity for small angles of propagation (un- 
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shaded area), corresponding to large angles of incidence, and it is close to minus unity for 
large angles of propagation (dark area), corresponding to large angles of incidence. This 
implies that there is no significant reduction in polarization on reflection for reflection at 
sufficiently small angles, and that there is almost complete reversal of the polarization on 
reflection at sufficiently large angles. Strong depolarization occurs for reflection near the 
contour Pr = 0. There are two sections of the contour pr = in figure 1, one in the region 
of total internal reflection and the other in the region where both reflection and transmis- 
sion occur. Only the case of total internal reflection is likely to be important in practice, 
and then, for semiquantitative purposes, this contour may be approximated by (10), which 
corresponds roughly to a straight line from uj = uj p , 9 = to uo = uj' p , 9 = 45°. Any reflection 
that occurs near this line leads to almost complete depolarization. 

The region where total internal reflection occurs is separated from the region where 
both reflection and transmission occur; in the two contour plots in figure 1 these regions are 
to the lower left and the upper right, respectively, of the curve corresponding to pr — 1. 
The region below this curve, where total internal reflection occurs, is restricted to smaller 
and smaller angles, 9, as uj/uj p increases. Although such small-angle reflection might be 
effective in causing ducting, it is not effective in causing depolarization. Above the curve 
Pr = 1 in figure 1, both reflection and transmission occur at the boundary. For reflection 
to be effective in causing depolarization in this case two conditions need to be satisfied: the 
depolarization coefficient should be close to zero, and the reflection coefficient, r tot , should 
not be too small. The depolarization coefficient well to the right of the curve p R — 1 in 
figure 1 is well approximated by setting 9' = 9 in (6), which reproduces the formula given 
by Wentzel, Zlobec and Messerotti (1986). In this regime, the contour line for Pr = 
asymptotes to 9 — > 45° for uj/uj p — > oo, corresponding to a deflection through 29 = 90°. As 
argued by Wentzel, Zlobec and Messerotti (1986), it is difficult to account for the observed 
depolarization in terms of a single reflection with this depolarization coefficient, and multiple 
reflections are required to account for substantial depolarization. Each reflection decreases 
the reflected intensity by a factor equal to r tot , and if r tot is very small multiple reflections 
lead to an extremely small intensity that can be of no interest in practice. The reflection 
coefficient is plotted in figure 2. One has r tot = 1 for total internal reflection, and figure 2 
shows that r tot falls off very rapidly in both frequency and angle above the threshold where 
transmission is allowed. Thus, in the regime well to the right of the curve pr = 1 in figure 1 
the intensity of the reflected component is orders of magnitude smaller than the transmitted 
intensity. It is concluded that effective depolarization requires total internal reflection, and 
that only the region to the left of the curve pr — 1 in figure 1 is likely to be relevant in 
developing a model for depolarization. 
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3.5. Validity of the model 

The model used here to calculate the depolarization coefficients corresponds to the 
limiting case of the magnetoionic theory in which the ratio of the electron cyclotron frequency 
to the plasma frequency is assumed negligibly small. The actual value of this ratio is thought 
to be small, of order Vt e /uj p pa 0.1 in type III sources, but may be larger in type I sources. 
Thus the formulae derived here are expected to be accurate only to of order 10%, which 
should suffice for semiquantitative purposes. However, there are special cases where the 
approximation fails and might be misleading, specifically for uj — uo p < Q e and uj — uo' p < Q e . 
In these cases, the refractive indices and the polarization vectors of the magnetoionic modes 
depend sensitively on the ratio of the frequency differences, uj — uj p , uj — u' p , to Q e . In 
particular, there is a range of frequencies above uj p (uj' p ) on the lower (higher) density side 
where the o mode can propagate and the x mode is evanescent. There is also a range around 
of frequencies ~ Q e around uj' p in the denser medium where a transmitted component in the 
z mode can propagate. 

The invalidity of the model for uj ~ uj p implies it cannot be used to treat any reflection 
that occurs near the point of emission. The x mode exists only at frequencies above its 
cutoff, uj > uj p + f2 e /2, so that there can be no reflected x mode for uj /uj p < Q c /2uj p . If one 
defines the reflection coefficients r ox , r OG for incident o mode and reflected o and x modes, 
respectively, then one has r ox = below the x mode cutoff. Moreover, there is a range of 
slightly higher frequencies where the refractive index of the x mode is much smaller than that 
of the o mode, and this mismatch implies r ox <C r OQ . Thus the approximation in which the 
effect of the magnetic field on the refractive index is neglected is valid only for uj — uo p ^> Q e . 
This limitation is unimportant provided the first reflection occurs well above the source 
region, where both the o and x modes can propagate, and this is assumed to be the case in 
the following discussion. 

The accuracy of the approximate model assumed here can be tested only by comparison 
with calculation using the exact magnetoionic expressions. Numerical calculations in the 
general case have been presented by Hayes (1985a,b). An example of the exact case is 
illustrated in figure 3. This case is similar to that shown on the left in figure 1 for 9 pa 60°. 
The exact result shown in figure 3 is for an incident x mode, and the results show that 
the polarization approximately reverses sense over most of the frequency range below the 
threshold for transmission: specifically, the reflected o mode is much stronger than the 
reflected x mode. However, this is not the case near uj = ui p and near the threshold for the 
appearance of a transmitted component. As already remarked, the approximate theory is 
invalid near uo = uj p , and it fails to predict the correct behavior near uj = uo p . However, 
it does predict the sharp narrow peak with pa 1 near the threshold for transmission. 
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This suggests that the approximate model is valid for semiquantitative purposes except near 

UJ = Up. 

The calculations performed here are all for the case where the magnetic field, the incident 
ray and the normal to the boundary are coplanar. This limitation is unimportant and 
the formulae apply to the more general case, provided that the angle 9 is reinterpreted 
appropriately. It is relatively simple to show this in the high-frequency limit (Appendix B); 
for arbitrary angles 9, 0, the result implies by (6), with 9' = 9, is reproduced with 29 
interpreted as the angle of deflection of the ray in this more general case. The foregoing 
results are insensitive to the assumption of coplanarity, and their generalizations involve 
only a reinterpretation of the angle 9. 

4. Depolarization of F emission 

In applying the foregoing results to the interpretation of F emission, it is important to 
distinguish between type I emission and type III emission, which have different characteris- 
tics. (The polarization of F type II emission is not discussed separately here: it is similar 
to F type III emission.) Type I emission from sources near the central meridian experience 
little depolarization, and the degree of depolarization increases systematically with displace- 
ment of the source towards the limb. Type III is never 100% polarized, with a maximum 
polarization of about 70%, with some bursts having a much lower polarization. A working 
hypothesis is that these differences are due in part to different density structures. Type I 
emission is assumed to come from closed-loop regions where the magnetic field is relatively 
strong and the density structures are aligned along field lines whose direction changes only 
over distances of order the major radius to the loop. Type III emission is assumed to come 
from open field regions, where the magnetic field is weaker and the density structures are 
aligned along field lines that have a random component about a mean orientation. 

4.1. Effect of refraction into the radial direction 

An important effect for any F emission is an initial strong refraction towards the local 
direction of the outward density gradient, which for the purpose of discussion is assumed to 
be the radial direction. The theory of F emission implies that the radiation is generated at 
a frequency that is above the local plasma frequency by fraction f ps 3V e 2 /2f^, where V e is 
the thermal speed of electrons and is the phase speed of Langmuir waves, identified as 
the speed of the electrons exciting the emission. For flare-associated type III emission this 
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fraction is of order 0.01, and is somewhat larger for storm associated type III emission, and 
perhaps also for type I emission. The refractive index at the point of emission is n V^fo- 
Let ip be the angle between the ray and the radial direction. Refraction implies that n sin if) is 
constant along a ray, Hence, along a ray, if) varies so that it satisfies (1 — uo 2 /u 2 ) sin 2 if> = 2/q . 
For example, for /o = 10~ 2 all the rays are confined to within ip 10° of the radial direction 
when the plasma frequency has decreased by a few percent. For magnetic field lines roughly 
in the radial direction, it then follows that when a ray encounters a boundary the angle 9 is 
determined primarily by the angle between the magnetic field and the radial direction. 

Consider the implication in connection with the depolarization coefficients shown in 
figure 1. F emission is generated very near the left-hand boundary of the figure. Radial 
refraction causes it to move rapidly to smaller angles as the frequency increases. If the ray 
does not encounter a boundary, it moves further towards the lower right of the figure, and if 
there were no other effect, the ray would emerge within a few degrees of the radial direction. 
Either refraction or reflection off density irregularities must occur to produce the inferred 
relative broad angular patterns of escaping radiation. 

4.2. Type I emission 

The model of Wentzel, Zlobec and Messerotti (1986) accounts for most of the important 
features of the directivity and depolarization of type I sources as they approach the limb. In 
the model, before encountering the boundary, the radiation is assumed to fill a cone of angles 
about the radial direction, and the magnetic field line that defines the boundary is at an 
angle 9 to the radial. Wentzel, Zlobec and Messerotti (1986) argued that the depolarization 
must take place well above the source region, specifically where uj/oo p > 2, in order to be 
consistent with a constraint on the time delay between the o-mode and x-mode components. 
This model is adopted here, except that the reflection is attributed to the density jump 
at the surface of an overdense fiber, rather than to scattering by lower-hybrid waves. The 
important change is in the depolarization coefficient, which has different forms for these two 
types of reflection. 

The depolarization coefficient for scattering by lower-hybrid waves corresponds to (6) 
with 9' — > 9. An attractive feature of this model is that the depolarization depends only on 
the geometry, specifically, the angle through which the ray is deflected. It correctly predicts 
unpolarized radiation for 29 = 90°, corresponding to a source on the limb. However, the 
model encounters a difficulty: averaging the reflected radiation over the assumed cone of 
angles makes it difficult to account for unpolarized type I sources near the limb. This led 
Wentzel, Zlobec and Messerotti (1986) to assume multiple (at least three) reflections in order 
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to account for the depolarization. For reflection at a sharp boundary, the same formula for 
the depolarization applies for uj ^> uj' p . However, the reflection coefficient, r tot , is then 
very small, cf. figure 2, and after three reflections the intensity would be negligibly small. 
To account for effective depolarization without a large reduction in intensity, the reflection 
must occur near or below the threshold for transmission to be allowed. The only plausible 
possibility is thus depolarization on total internal reflection. 

Three constraints severely limit the possibilities: the need for total internal reflection, 
the inference from the time delay between the o and x mode component implying depo- 
larization at uj/uj p > 2, and the requirement that the deflection occur through roughly 90° 
(0 « 45°) for sources near the limb. These constraints can be satisfied simultaneously only 
for a relatively large density ratio, £ ~ 10, with the depolarization coefficient is then given 
by (9). It is possible to have a relatively large depolarization for a single reflection. How- 
ever, as in the model of Wentzel, Zlobec and Messerotti (1986), for radiation in a cone of 
angles it is difficult to account for very small degrees of polarization with a single reflection. 
There is an additional effect that is absent in the model of Wentzel, Zlobec and Messerotti 
(1986): the depolarization depends on frequency as well as angle. According to figure 1, a 
given fiber reflecting radiation with a range of frequencies at a given angle has pn increasing 
with frequency; if Pr = occurs in this frequency range, then the lower frequency radiation 
should have the opposite handedness to that of the higher frequency radiation. Observation 
of such an effect would support the form of the model suggested here. 



4.3. Type III emission 

In contrast to type I burst, which are typically 100% polarized, F emission in type III 
bursts is never 100% polarized. In the framework of the present model, this requires that 
the radiation always experience at least one depolarizing reflection before escaping. 

The picture of the source region assumed here is a low-density background region, with 
overdense field-aligned structures acting as the walls of ducts (Duncan 1979). For simplicity 
let us assume that all these structure have the same overdensity, £. The structures are also 
likely to have a distribution in angle, 9, due to random wandering of the open field lines. 
The near-radial ducting suggests that the mean value of 6 is relatively small: in the following 
discussion a mean 9 « 20° is assumed. Type III emission is assumed to be generated in the 
low-density region, to be refracted quickly towards the radial direction, to experience its first 
reflection off one of the overdense structures, to experience at least a few more reflections 
as it propagates outward, and to escape when it encounters no more structures capable of 
causing total internal reflection. The data suggest that both F and H emissions are ducted 
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(Duncan 1979), implying at escape occurs above the second harmonic plasma level, with 
u/ujp > 2 at the escape point. The condition for total internal reflection to cease at the 
escape point is 



Assuming uj/uj p ps 2.5 and 9 ps 20°, (11) suggests that structures with £ ps 2 are involved in 
ducting. 

Of particular importance for depolarization is the first reflection. According to figure 1 
for £ = 2, for 9 ps 20° all reflections at cj/c<jp > 1.5 have > 0.8, implying that many such 
reflections would be required to reduce the polarization even to the highest observed values, 
~ 70%. A possible model for depolarization of type III bursts is that the first reflection 
occur at 1 <C uo/uo p < 1.5, and leads to a significant depolarization, and that two or more 
subsequent reflections occur, leading to further depolarization by lesser amounts. For the 
most highly polarized bursts, the model suggests u/u p ~ 1.5 with a further reflection at 
u/ujp ~ 2 before escaping at uj/oj p ~ 2.5. Lower polarization is more plausibly attributed 
to the first reflection being closer to the source, leading to a greater depolarization, rather 
than to a greater number of subsequent reflections, which become increasing ineffective in 
causing depolarization. 

It may be concluded that depolarization due to reflection off the walls of a duct can 
account for the observed polarization of F type III bursts. The walls of an appropriate 
duct would consist of overdense structures with £ ps 2 along magnetic field lines oriented at 
9 pa 20° to the radial. The model suggests that all type III bursts experience at least one 
reflection, with several further reflections likely before escape from the duct at uo/uo p > 2. 
Most of the depolarization occurs at the first reflection. 



Effective depolarization requires that radiation initially in one magnetoionic mode be- 
comes a mixture of the two magnetoionic modes, and although reflection is one possibility, 
there is another possibility. Mode coupling is a weak effect except under two conditions: 
when the density gradients are very large, as assumed here, and near coupling points. In 
the absence of sharp density gradients, mode coupling in the corona is important only at 
so-called QT regions, which are coupling points where the angle, 9, between the ray and 
the magnetic field passes through 90°. (QT is an abbreviation for 'quasi-transverse' in now 
outdated jargon: in modern usage this corresponds to 'quasi-perpendicular' when the po- 
larization is better approximated by linear rather than circular.) When a QT region is 
encountered at sufficiently low frequencies, the radiation remains in the initial magnetoionic 




(11) 



4.4. Depolarization at a QT region 
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mode and the sense of polarization reverses, whereas for an encounter at sufficiently high 
frequencies, the sense of polarization is preserved, (Budden 1961; Zheleznyakov 1970; Melrose 
1980). Significant depolarization at a QT region occurs only at frequencies near the tran- 
sition frequency between these limits (e.g. Melrose 1980, vol. II, p. 300). In principle, this 
provides an alternative depolarizing mechanism. 

No model for depolarization of F emission has been proposed based on QT regions. Any 
such model would encounter two major difficulties, based on the requirements that a QT 
region be encountered near the transition frequency. The geometry required for the rays 
to encounter an appropriate QT region does not fit with either the fiber model for type I 
sources or the ducting model for type III sources. Furthermore, there is no obvious reason 
why in any specific model a QT would be encountered near the transition frequency. 

5. Conclusions 

In this paper it is assumed that depolarization of F emission in type I, II and III bursts is 
due to reflection. In general, a wave in one magnetoionic mode incident on a sharp boundary 
leads to reflected waves in both modes, and to transmitted waves in both modes provided 
these waves can propagate in the denser (by a factor £) medium. The reflection model 
requires that the boundary have a thickness less than about the (free-space) wavelength of 
the radio waves. 

The depolarization is treated in section 3 in terms of a depolarization coefficient, pr, 
with p R = corresponding to complete depolarization, and Pr — ±1 corresponding to no 
change in polarization and reversal of (circular) polarization, respectively. A simplified form 
of magnetoionic theory is used to calculate pr. Two different expressions are derived for 
Pr. (6) applies when the transmitted waves can propagate, with the reflection coefficient 
then given by (7), and (9) applies when total internal reflection occurs. Contour plots of 
Pr as a function of uo/uo p and 9 are shown in figures 1: the general form is not sensitive 
to £, provided the frequency is scaled appropriately, as exhibited by the approximate form 
(10). The simplified form of magnetoionic theory reproduces relevant features found from 
exact calculations of reflection coefficients (Hayes 1985a,b), except for uo u p ,lu' p , where 
the neglect of the magnetic field on the refractive indices of the magnetoionic waves is not 
justified. 

In applying the model to F emission, it is assumed that the radial gradient in the corona 
density causes the ray to be strongly refracted towards the outward radial direction before 
it encounters a sharp boundary. The invalidity of the simplified model for uj « uj p is then 
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unimportant, and the angle 9 may be interpreted as approximately the angle between the 
magnetic field and the radial direction. It is found that when transmission is allowed the total 
reflection coefficient is very small except for a narrow range of parameters that is unlikely to 
be important in applications. On this basis, it is argued that effective depolarization requires 
total internal reflection. 

The depolarization of type I is attributed here to reflection off the overdense structures 
(fibers) invoked to account for the directivity (Bougeret and Steinberg 1977). Type I emission 
from sources near the central meridian must escape without encountering such a structure. 
For sources nearer the limb, the model proposed by Wentzel, Zlobec and Messerotti (1986) 
is adopted, except that the reflection (that these authors attributed to scattering off lower- 
hybrid waves) is assumed to be off a sharp boundary. The time delay between the o and 
x mode components requires reflection well about the source region (Wentzel, Zlobec and 
Messerotti 1986), and for this to correspond to total internal reflection through a large angle 
(tending to 29 = 90° for a source at the limb) requires a large density contrast, £ ~ 10. 
Depolarization of F emission in type II and III bursts can be explained in terms of reflection 
off ducts with a density contrast £ ~ 2. A ray then typically experiences two or more 
reflections before escaping, with the first reflection at uj/uj p m 1.5. 

An important implication of the reflection model for depolarization is that there must 
be density structures with much more extreme properties than required by observations at 
non-radio wavelengths. There are few alternatives explanations for depolarization that would 
avoid this implication. One alternative is that proposed by Wentzel, Zlobec and Messerotti 
(1986) for type I bursts: scattering off lower- hybrid waves generated through some appro- 
priate instability. The lower-hybrid model requires that the lower-hybrid waves be present 
along a thin boundary layer that extends over a macroscopic distance and persists for many 
hours, and this is difficult to reconcile with the typical burstiness, in both space and time, 
of instabilities. Scattering by lower-hybrid waves has not been suggested for depolarization 
of F emission in type II and type III bursts, and such a suggestion would encounter serious 
difficulties. In particular, the angles of deflection are thought to be relatively small dur- 
ing ducting, leading to a correspondingly small depolarization, whereas the depolarization 
is universal and relatively large for these bursts. Another possible depolarizing mechanism 
is associated with mode coupling at so-called QT regions. This mechanism has not been 
proposed for depolarization F emission, and any model based on it would require very spe- 
cial conditions to apply. Another possibility that might be considered is that depolarization 
occurs very near to the source due to some unknown process that coverts o mode waves into 
x mode waves. Even if one postulated such a process, it would be inconsistent with the 
measured time delay between the o and x mode components in partially polarized type I 
source, which imply depolarization at u/uj p > 2 (Wentzel, Zlobec and Messerotti 1986). Thus, 
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the only seemingly viable depolarization mechanism for F emission is reflection off density 
structures with large density contrasts and extremely sharp boundaries. 

I thank Qinghuan Luo and Mike Wheatland for comments on the manuscript. 

A. Condition for reflection 

Consider a ray encountering a boundary layer in which the plasma parameters vary as a 
function of only the coordinate z. Snell's law implies that r = (k 2 x + ky)^ 2 c/u is a constant, 
and it is convenient to solve the magnetoionic dispersion equation (the 'Booker quartic' 
equation) for the dependent variable to q — k z c/uj. A ray propagating from the lower 
density towards the higher density region encounters the reflection point for the o mode 
before the reflection point for the x mode. The condition for reflection to apply, in the sense 
that a continuous boundary layer is treated as a discontinuity, is that the skin depth for 
the o mode in the region of evanescence between these reflection points exceed the distance 
between the reflection points. For Q e <C u p , the frequency u> is approximately equal to the 
plasma frequency, uj p near the reflection points, which occur at 1 - X - Y - r 2 w and 
1 — X — r 2 ps 0, respectively. If L is the characteristic length over which the density changes, 
the reflection points are separated by Az ps (Vl e /uo p )L. The skin depth is identified as the 
inverse of (u/c)lmq. The relevant approximate solution of the Booker quartic equation (e.g. 
Budden 1961) is for Q e /u <C r, and it may be approximated by 

Imq = (n e /w p )[(l -r 2 )(l + r 2 cos 2 0)] 1/2 , (A.l) 

where tp is the angle between the magnetic field and the z-axis. Assuming the magnetic 
field lines are parallel to the surfaces of constant density implies cos^ = 0, and then the 
condition for reflection becomes (1). 

B. Depolarization at high frequencies 

The following derivation reproduces the reflection coefficient written down by Wentzel, 
Zlobec and Messerotti (1986), and indicates how to relax the coplanarity assumption made 
in section 3. 

Consider a model in which the magnetic field is along the 3 axis and the normal to 
the boundary is along the 1 direction. Let the incident wave be at polar angles 0, rel- 
ative to the magnetic field, with = corresponding to the 1-3 plane. The incident 
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wave is then along the direction (sin 9 cos 0, sin 9 sin 0, cos 9). The reflected wave is along 
(— sin# cos0, sin 6* sin 0, cos 9). Right and left hand polarizations for the incident wave are 

e^L = —i=(cos 9 cos =|= 2 sin 0, cos 9 sin ± 2 cos0, — sin#). (A.l) 
v2 

Circular polarizations for the reflected waves, e'^ L say, are given by cos0 — > — cos0 in (A.l). 
The reflection coefficients are in the ratio 

trr : f'RL = \e* R ■ e' R \ 2 : \e* R • e" L \\ (A.2) 

where an asterisk denotes complex conjugation. A straightforward calculation then gives 

(cos 2 9 cos 2 + sin 2 0) 2 — (sin 2 9 cos 2 0) 2 
PR ~ (cos 2 9 cos 2 + sin 2 0) 2 + (sin 2 9 cos 2 0) 2 ' ^ ' ' 

The result (A. 3) simplifies in two special cases. For radiation incident in the 1-3 plane 
(0 = 0) one has p R = (cos 4 # — sin 4 9) /(cos 4 9 + sin 4 #), and for radiation incident in the 
1-2 plane {9 = vr/2) one has pn = (sin 4 — cos 4 0)/ (sin 4 + cos 4 0). These two result are 
identical when written in terms of the angle of incidence (with the respect to the normal to 
the boundary), 9 inc , which is n/2 — 9 and 0, respectively. Inspection of (A. 3) shows that the 
general result may also be expressed in this same form, that is 



sin 4 9 mr — cos 4 9\ nc 2 cos 9. 



" jiu ' ' ' "Die • •* (irf /a , \ 

Sin fc'inc + COS 4 6> inc 1 + COS fldef 

where ^dcf = vr — 29 inc is the angle through which the ray is deflected. The result (A. 4) was 
written down using a more qualitative argument by Wentzel, Zlobec and Messerotti (1986). 

The derivation in section 3 of the reflection coefficient off a sharp boundary is for = 0. 
The generalization to ^ for simply by reinterpreting 29 as #def- 
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Fig. 1. — The depolarization coefficient is shown as a contour plot on the uj/uo p -9 plane 
for £ = ujp/ujp = 10 (left) and 2 (right). The contours are, from bottom to top, p R = 
0.8, 0.4, 0, —0.4, —0.8. Thus the unshaded region correspond to depolarization coefficients 
0.8 < pr < 1, and the dark region corresponds to — 1 < p R < —0.8. The curve contained 
within the unshaded region, and between the contours where they converge on the upper 
right, corresponds to pr — 1; this curve separates total internal reflection, to its left, and 
partial reflection and partial transmission, to its right. The ranges of frequency for the two 
figures are chosen to emphasize the approximate scaling as a function of £ 1//2 . 




Fig. 2. — The reflection 
function of u/u p for 6 = 



coefficient r tot is plotted for a density ratio £ = 10 is plotted as a 
45° (left) and as a function of 9 (in degrees) for u/u p = 5 (right). 
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Fig. 3. — The reflection coefficients (S z is the reflected Poynting flux along the z axis for 
unity initial Poynting flux equal to unity) and the solution, q 2 of the Booker quartic equation 
are shown for an incident x mode on a boundary with a density contrast uo 12 /uo 2 = 8 with 
Vt e /uo p = 0.2 and 9 = 60°. The reflected x and o modes are shown on the left and right, 
respectively. [From Fig. 2 of Hayes (1985a).] 



